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Abstract 

A quantum system at equilibrium is represented by a corresponding classical system, chosen 
to reproduce thermodynamic and structural properties. The motivation is to allow application 
of classical strong coupling theories and molecular dynamics simulation to quantum systems at 
strong coupling. The correspondence is made at the level of the grand canonical ensembles for the 
two systems. An effective temperature, local chemical potential, and pair potential are introduced 
to define the corresponding classical system. These are determined formally by requiring the 
equivalence of the grand potentials and their functional derivatives. Practical inversions of these 
formal definitions are indicated via the integral equations for densities and pair correlation functions 
of classical liquid theory Application to the ideal Fermi gas is demonstrated, and the weak coupling 
form for the pair potential is given. In a companion paper two applications are described: the 
thermodynamics and structure of uniform jellium over a range of temperatures and densities, and 
the shell structure of harmonically bound charges. 



I. INTRODUCTION AND MOTIVATION 



The equilibrium properties (thermodynamics and structure) of classical systems for condi- 
tions of strong coupling (e.g., liquids) can be described by a number of practical and accurate 
approximations pp. Examples are those based on integral equations for the pair correlation 
function, such as the Percus-Yevick and hypernetted chain approximations (HNC). Ther- 
modynamic properties are then determined from exact representations in terms of the pair 
correlation function. Extensions to inhomogeneous systems have been described in terms of 
the pair correlation function as well. More extensive formulations of the thermodynamics 
and structure for classical systems can be obtained within the framework of classical density 
functional theory [2] . Molecular dynamic simulation (MD) is perhaps the most accurate tool 
for the description of classical systems. 

A corresponding practical description of strong coupling when quantum effects (diffrac- 
tion and degeneracy) are important is less well-developed, except when a transformation to 
a corresponding weak coupling quasi-particle representation can be found (e.g., phonons, 
conventional Fermi liquid theory). Conditions of current interest for which such representa- 
tions are not available include those of the growing class known as "warm, dense matter", 
i.e. complex ion-electron systems [3]. The ions are typically semi-classical but the electron 
conditions span the full range of classical plasmas to zero temperature solids. With such 
potential applications in mind, the objective here is to extend the classical approximation 
methods to include quantum effects. The idea is to introduce a representative classical 
statistical mechanics that embodies selected quantum effects via effective thermodynamic 
parameters and effective pair potentials in the classical Hamiltonian. These classical pa- 
rameters and pair potential are defined formally in such a way as to ensure the equivalence 
of the classical and quantum thermodynamics and structure. In this way the established 
effective classical methods can be applied to describe quantum systems at strong coupling 
as well. The precise definition of the effective classical statistical mechanics is given in the 
next section. Of course, there is no possibility to map a quantum system entirely onto an 
equivalent classical system. Instead, the classical system is defined to yield only selective 
quantum properties, in this case the thermodynamics and pair structure. Other properties 
calculated from this effective statistical mechanics, such as higher order structure or trans- 
port properties, constitute uncontrolled approximations, and its utility in this more general 
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context must be obtained from experience with applications. 

The analysis here is a formalization of related partial attempts to incorporate quantum 
effects in classical calculations, most commonly through effective classical pair potentials. 
There is an extensive history for constructing such effective potentials [IH9] . Weak exchange 
symmetry (degeneracy) effects were first described via a classical pair potential by Uhlenbeck 
and Gropper [3], and were extended to arbitrary degeneracy for the ideal (non-interacting) 
quantum gas by Lado more than fifty years ago [7j. Pair potentials representing diffraction 
effects have typically been defined in terms of the two particle density matrix [S]. The 
combined degeneracy and diffraction effects are typically combined additively into a single 
pair potential pU]. Such potentials have been applied in a number of classical theories and 
simulations. In particular they provide a definition for the classical statistical mechanics 
of electron-ion systems which otherwise do not exist due to the Coulomb collapse without 
diffraction effects. A related approach resulting in a classical dynamics with an effective pair 
potential is "wave packet molecular dynamics" [T2] . This variational approach is not based 
in thermodynamics and is somewhat outside the topic of discussion here. 

Effective pair potentials defined from the two particle density matrix do not assure that 
the thermodynamics or structure will be given correctly by the corresponding classical sta- 
tistical mechanics, and those approaches provide uncontrolled approximations in most ap- 
plications. A more general method defines a classical Boltzmann-Gibbs ensemble from the 
N particle density matrix, or Slater sum [HI HH H2]- This gives not only the correct ther- 
modynamics and pair structure, but all higher order multi-point correlations as well. The 
price of such a definition is that the effective classical potential is no longer given in terms of 
pair potentials, but includes all higher order many particle potentials as well. Applications 
typically truncate this at the level of the pair potentials, leading again to an uncontrolled 
approximation. An equivalent Feynman-Kac functional representation has only pair inter- 
actions but among additional "internal degrees of freedom" (point particles are replaced 
by extended filaments or polymers whose "beads" interact) [Hj. In contrast, the formal 
approach here has the advantage of predicting the correct thermodynamics and structure 
while retaining the simplicity of pair potentials for point particles without approximation. 
However, the additional complication implied by this is the need for effective classical ther- 
modynamic parameters (e.g., temperature and chemical potential) that differ from those for 
the given quantum system. 
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The introduction of an effective temperature in this context appears to have been done 
first by Perrot and Dharma-wardana (PDW) in their classical map hypernetted chain model 
[IB*] . It provides a means to calculate a pair correlation function g c (r) from the classical HNC 
integral equation with two modifications: 1) the actual pair potential is replaced by the sum 
of a pair potential to give the correct ideal gas quantum pair correlation function, plus a 
phenomenological modification of the actual pair interaction to account for diffraction effects 
(Deutsch potential [IB]), and 2) an effective temperature T c = (T 2 + TJ 2 ) 1//2 . The single new 
parameter T is determined as a function of the density by requiring that the classical 
correlation energy at T = be the same as that for the quantum system. The latter is 
determined from diffusion Monte Carlo simulation data. The solution to the HNC integral 
equation with these embedded quantum effects combined with classical strong coupling 
provides an approximate g c {r) from which thermodynamic properties can be calculated from 
their classical expressions (e.g., pressure from the virial equation). During the past decade 
this remarkably simple and practical approach has been applied to a number of quantum 
systems with surprising success for a wide range of temperatures and densities. A recent 
review with references is given in [17J. These successes provide motivation and support for 
the work presented here, which can be considered a formalization of these concepts into a 
precisely defined context. 

The definition of a representative classical system in terms of selected equilibrium prop- 
erties of the underlying quantum system is given in the next section. The quantum system 
is represented in the grand canonical ensemble. For generality an external single particle 
potential 4> ex t (r) is included so the thermodynamics is that for an inhomogeneous system. 
Interactions among particles is given by a pair potential <fi (r, r'). Only a single component 
system is considered here for simplicity. Thermodynamic properties are functions of the 
temperature T and functionals of the local chemical potential /x (r) = /x — <p ext (r). A corre- 
sponding classical system is specified in a classical grand canonical ensemble in terms of an 
effective classical temperature T c , classical local chemical potential /x c (r) = /x c — (p c ^xt (r), 
and classical pair potential c (r,r'). These three classical quantities are defined by three 
independent constraints equating chosen classical and quantum thermodynamics and struc- 
ture: equality of the grand potentials and their functional derivatives with respect to the 
chemical potential and pair potential for the classical and quantum systems. Equivalently, 
this requires that the the pressures, the local densities, and the pair correlation functions 
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are the same for the two systems. 

Having defined the classical system in terms of equating classical and quantum proper- 
ties, there is then the practical problem of inverting the defining equations for the classical 
temperature, local chemical potential, and pair potential. Formally this is done via classical 
density functional theory which relates the local density and pair correlation function to the 
classical potential. Inverting this and exploiting the equivalence of classical and quantum 
pair correlation functions and densities gives the desired expression for the classical pair 
potential as a functional of the quantum pair correlation function and density. Examples 
of practical implementations of this approach are given by Percus-Yevick (PY) and HNC 
approximate integral equations. With the classical pair potential determined (formally), the 
classical temperature is obtained from the classical virial equation for the pressure and the 
HNC equation for the chemical potential. 

In this way the classical temperature, chemical potential, and pair potential are given 
as explicit functionals of the quantum properties. In practice, analysis of these functionals 
entails the full quantum many-body problem so it would seem that little progress has been 
made. The key assumption in this approach is that the strong coupling effects have a 
dominant classical component while dominant quantum effects are local (e.g. diffraction) 
or global but weakly dependent on interactions (e.g. exchange degeneracy). Then, simple 
weak coupling approximations to the defining functionals for the classical parameters are 
sufficient while strong coupling effects are obtained by a more complete implementation of 
the classical statistical mechanics. A more ambitious approach would try to include limited 
strong coupling effects in the classical parameters (e.g. T = correlation energy as in the 
PDW model above). The formalism here provides a basis for approximations that can be 
tailored to the specific system or objective considered. A brief overview of this approach has 
been described elsewhere [18] 

As a first application the uniform electron gas (jellium) is considered in the following 
companion paper [19] . The effective pair potential used incorporates the exact ideal gas 
exchange and the random phase weak coupling limit. The pair correlation function is then 
obtained for this potential from the HNC integral equation describing strong classical corre- 
lations, known to be accurate for Coulomb interactions. The results are compared to earlier 
quantum methods |20j and to the above PDW model over a range of temperatures and 
densities. Good agreement with diffusion MC simulations at T = is obtained. A second 
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application described briefly there is to charges in a harmonic trap where classical strong 
correlations produce shell structure. Qualitative effects of diffraction and exchange on the 
formation of shell structure is observed. 



II. DEFINITION OF THE REPRESENTATIVE CLASSICAL SYSTEM 
A. Thermodynamics from statistical mechanics 

Consider a quantum many-body system with Hamiltonian if in a volume V at equilibrium 
described by the grand canonical ensemble, with inverse temperature (3 and chemical po- 
tential \i . For simplicity a one component system is considered and dependence on internal 
degrees of freedom such as spin is suppressed. The Hamiltonian H is of the form 

H = K + $ + $ ext , (1) 

where K denotes the total kinetic energy, $ is the pair-wise additive potential energy among 
particles, § ext is an external potential coupling to each particle 

j N N 

$=2^^'^' ®ext = ^2<f>ext(<li), (2) 

ij i=l 

and N is the particle number. The statistical density operator for the grand ensemble 
depends on H in the form H — fiN. It is convenient to combine the contributions from the 
external potential with that from N 

H-iiN = K + $-J dr^(r)n(r), (3) 

where n(r) is the number density operator 

N 

n(r) = (4) 

i=i 

and /i(r) is a local chemical potential 

M r ) = A* - ^ext (r) . (5) 

A caret over n is used to distinguish the operator from its average value n. 

The thermodynamics of the system is determined from the grand potential Q = —f3~ 1 In S, 
where S is the grand partition function. The thermodynamic parameters for this function 
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are the inverse temperature j3, the local chemical potential /x(r), and the volume V. The 
grand potential is 

Q(f3 \n,4>) = -/3- 1 \nE(f3 | //, 0) = -/T 1 In ^ Tr* e^(* + *-^ (r) " (p) ) . (6) 

N 

The symbol Tr N denotes a trace over the N particle Hilbert space with appropriate sym- 
metry restrictions for Fermions or Bosons. The vertical line in the arguments of Q and S 
denotes a functional of the quantity following it, e.g. S(/3 | /i, 0) is a function of (3 and a 
functional of /x(r) and 0(q;, q^). Its functional dependence on the pair potential 0(r, r') has 
also been made explicit in this notation, although it is not strictly a thermodynamic vari- 
able. The dependence on the volume V is left implicit. The pressure p((3 \ ji) is proportional 
to the grand potential 

vip\p y <l>)v = -n{p\p,4>). (7) 

The first order derivatives of the grand potential give the internal energy E(/3 \ fi, 0) and 
average number density n(r;/3 | //, 0) 

W^mzUhV^, (8 ) 

nW|M)V = -^|^|,. (9) 

Higher order derivatives provide the fluctuations (susceptibilities) and structure. In partic- 
ular, the second functional derivative with respect to /x(r) is related to the response function 
x(r,r';/3 | /i) 

1 5 2 Q{(3 | /v 



(3 5//(r)5//(r') 



1/3= -x(r,r;/3 | fx, 
1 /* 



J dp' ^' H 5n(r)e~^' H n(r'y, f3 \ //, 0^ , (10) 



where <5n(r) = n(r) — n(r), and (X; /3 | /z, 0) denotes an equilibrium grand canonical average 
of the quantity X. A related quantity is the pair correlation function, g(r,r';/3 \ 0), 
obtained by functional differentiation with respect to the pair potential 

5Q(f3 | /i,0 



<50(r,r') 



p,n = n(r;p \ /i, 0)n(r';/3 | fi, 0)#(r, r';f3 \ (i, 0) 

= (n(r)n(r'); \ ^ 0) - n(r;/9 | /i, 0)5 (r - r') (11) 
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For the quantum system there is no simple relationship between the two measures of struc- 
ture, x(r, r 7 ;/5 | /i, 0) and g(r,r';/3 | Instead, only their time dependent extensions are 
related via a fluctuation - dissipation relation. 

A corresponding classical system is considered with Hamiltonian H c in the same volume V 
at equilibrium described by the classical grand canonical ensemble, with inverse temperature 
0c and local chemical potential /x c (r) . The Hamiltonian has the same form as ([I]) except 
that the potential energy functions ^ are different, and denoted by 

N N 
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(|> ' =iX^(qi,qj), ^c,ext = ^2^c,ext(qi) ■ (12) 



The local chemical potential is /i c (r) 

// c (r) = /i c - c>ert (r) . (13) 

The classical grand potential is defined in terms of these quantities by 

/3fi c (/3 c | /i c ,0 c ) = -mS c (/3c I A*e,&) = -ln]T — ^— / ^..dq^^-^^^). 

(14) 

1 /2 

Here, A c = (27r(3 c h 2 /m) is the thermal de Broglie wavelength associated with the clas- 
sical temperature. The integration for the partition function is taken over the N particle 
configuration space. 

The classical thermodynamics is determined in the same way as in ^ - Q 

p c (/3 c | /i c , <p c )V = -Q C (/3 C | He, 4>c)- (15) 

^c(/3 C | VcAc) = qq k^A* (I 6 ) 

n c (r;/3 c | /i c , C ) = — — \ Mc . (17) 

1 ^(/3 c |/i c ,0 c ) , _ , , „ , , 
K WJWr') = " Xc(r ' r A 1 ^'^ 

= -(a(r)n(r / );/3 c |/x c ,0 c ) c (18) 

5fi(/3 c | /i c ,( 



r,r' 



|/9 ,Ai e = n c{r](3 c I ii c ,<j>c)n c (r';/3c I <j>c)g c (r, r';/3 c | /i c ,0 c ) 

= (<5ra(r)n(r'); /3 C | /i c , C ) C - n c (r;/3 c | /i c , C )5 (r - r) . (19) 
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Here (X; (3 C \ fi c , C ) C denotes the corresponding classical equilibrium grand canonical av- 
erage of the quantity X. Note the additional constraint that the derivatives are taken at 
constant pair potential C . This is necessary because the determination of (3 C , /i c (r), and 
C for the quantum correspondence in the next subsection implies c (r, r') is a function of 
(3 C and a functional of /U c (r). This dependence is left implicit to simplify the notation. See 
Section below for further elaboration. 



B. Classical - quantum correspondence conditions 

The classical system has undefined ingredients: the effective inverse temperature, (3 C , the 
local chemical potential, /i c (r), and the pair potential for interaction among the particles, 
c (r, r'). A correspondence between the classical and quantum systems is defined by express- 
ing these quantities as functions or functionals of (3, /x(r), and 0(r, r'). This is accomplished 
by requiring the numerical equivalence of two independent thermodynamic properties and 
one structural property for the classical and quantum systems. The first two are chosen to 
be the equivalence of the grand potential and its first functional derivative with respect to 
the local chemical potential. 

ttdPc I He, 0c) = M(/3 I n, 0), JJTJr) S/4r) ^ ' ^ ' 

An equivalent form for these conditions can be given in terms of the pressure and density 
p c (/3 c | /i c , C ) = p(/3 \fi,4>), n c (r;/3 c \ /i c , C ) = n(r;/3 \ fi, 0). (21) 

These two relations provide two independent relations between of j3 c , /i c ( r ) and the physical 
variables (3 and /i(r). 

It remains to have a structural equivalence to relate the pair potential c (r, r') to 0(r, r'), 
which are two particle functions. This is accomplished by equating the functional derivatives 
of the grand potentials with respect to these pair functions 

5ttJA I /Xc, 0c) | _ 5£l(P | /i,0) 

50 c (r,r') l/3c ' Mc " 50(r,r') ^ " { ' 

More physically, this implies the equivalence of density fluctuations 

(^(r)n(r'); & I 0c) c = (Sn(r)n(r') ; (3 \ //, 0) . (23) 



9 



Finally, from the equivalence of densities and the definitions of pair correlation functions in 



(11) and (19) this third condition becomes 

9c(r, r'; (3 C \ \i c , <j> 



g{r,r';f3 \ /i, 



(24) 



In summary, the classical - quantum correspondence conditions are the equivalence of the 
pressures, densities, and pair correlation functions. 



It might seem that a natural alternative choice to (24) would be to equate the second 



functional derivatives of the grand potential with respect to the local chemical potential, or 
equivalently x c (r, r';/3 c | li c ,4> c ) = x( r , r ';/3 I f 1 ,^)- However, the classical response function 
has a singular contribution proportional to 5 (r — r') that is not present in the quantum 
response function. Hence their equivalence does not provide a simple mapping of the pa- 
rameters. In contrast, the classical and quantum forms for the pair correlation functions are 
similar and do not have this problem. 



In this way the three equations of (21) and (24) determine, formally, the classical param- 
eters (3 C , /i c , and (fi c (q) as functions of (3, and functionals of fi(r), and <fi (r, r') 



0c = PdP I At, 0), A* c = A*c(r;/3 | li, 



(r,r';/3 | li, 



(25) 



This completes the definition of the classical system representative of the given quantum 
system. 



C. Inversion of the correspondence conditions 



To determine the classical parameters, (25), inversion of the functional forms (21) and 



(24) is required. This is done on the basis of three exact results of classical statistical 



mechanics, for the pressure, density, and pair correlation function respectively. The first 
is the virial form for the classical pressure (obtained by differentiating Q c (f3 c | fi c , C ) with 
respect to the volume) 

PcPc{Pc I /i c , <fic) = ^ J dvn c (r) 1 + ■ V/3 c /i c (r) - ^ / dr'n c (r / )^ c (r, r')r' ■ V'/3 c c (r, r) 

(26) 

(Note that this is the average pressure over the entire system. A local pressure could be 
identified with the integrand of (26)). The second is an equation for the classical densities, 

In (n c (r) A^) = /3 c /i c (r) + J da j dr"cf\v, r" | an c )n c (r") . (27) 
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and the third is an equation for the pair correlation functions 

lng c (r, r') = -/3 c c (r, r') + J da J dv"c {2) (j, r" | n c + an c (g c - 1)) 

xn c (r")(g c (r",r')-l) (28) 

(2) 

The direct correlation function c c appearing in these equations is determined from g c via 
the Ornstein-Zernicke equations 



(& (r, r') - 1) = c c (r, r' | n c ) 

+ J dr"c c (r, r" | n c ) n c (r") (# c (r", r') - 1) . 



(29) 



The constant A c in (27) is the thermal de Broglie wavelength evaluated at the classical 
temperature, A c = (27ih 2 p c /m) 1 ^ 2 . See Appendix 



A 



for further elaboration. 



These equations provide the inversion in the following way. First, l\27h and (28) are solved 



for (3 c fi c (r) and /3 c c (r,r') as functionals of n c (r) and g c (r,r'). Then n c (r) and g c (r,r') are 



replaced by their quantum counterparts, n(r) and g(r, r'), according to the definitions (21) 

3d 

'A 



and (24), giving the desired exact expressions 

+ In (n (r) A 3 ) - J da j dr"c (2) (r, r" 



&/i c (r) 



?ln, 

2 V/3 



cm m r 



(30) 



/3 c c (r, r') = — In g(r, r') + J da J dr"c^ 2 \r, r" \ n + an (g — 1)) 

xn(r")(#(rV')-l). (31) 



Here, c^(r, r" | n) is defined by the Ornstein-Zernicke equation (29) with n c (r) , g c (r, r') 
replaced by n (r) , g(r, r') 



c ( 2 ) ( r r ' I n \ 



(g (r, r') — 1) — J dr"c^ (r, r" \ n) n (r") (g (r", r') — 1' 



(32) 



Finally, an equation for /3 C / /3 is obtained by using the equivalence of the classical and quan- 



tum pressures p c i = p, (21 ), to write 



Pd _ PclPcl 

P : Pp ' 



or with (26) 



1 



Pc_ 

P Ppv 



J dvn(r) 1 + ir • V/3 c/ u c (r) - ^ / dr'n(r')<?(r, r')r' • V'/3 c c (|r - r'|) 



(33) 



• (34) 
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Equations (30), ( |3T| ), (32), and (34) then determine the classical parameters as (3 c /f3, 
j3 c fj, c (r), and /3 c c (r,r') from the properties of the given quantum system 

P c = p c {P\ n, g) , /i c (r) =/i c (r,/3 | n,g), c (r, r') = c (r, r', p \ n, g) (35) 



which is equivalent to (J25| (corresponding to a change of variables from fi, <fi to n, g for the 
quantum system). 

These expressions for the classical parameters are exact, but deceptively simple, equa- 
tions. Their complexity arises from the fact that different arguments for the functional 
c^(r, r" | ■) occur. The relationship of two functionals at different arguments depends 
on higher order correlations that are not determined by these equations. The many-body 
problem has not been solved, only hidden, and approximations are required for practical 
applications 

The application of these results proceeds in three steps. First an approximation to the 



classical equations (27) and (28) is chosen such that classical correlations are well described. 



These are then inverted to obtain the corresponding approximations to (30) and (31 ). Next, 



selected information about the quantum system is chosen in the form of approximations to 



p, n, and g, for the calculation of the classical parameters of (35). Finally, with all classical 



parameters known properties of interest are determined from classical statistical mechanics 
(e.g., liquid state integral equations, classical density functional theory, molecular dynamics). 



D. Hypernetted Chain Approximation 



To illustrate the first step, selection of a classical approximation for (27) and (28), the 



HNC is noted. This results from making them local in the function space of densities 

cf\r, r" | n c + an c {g c - 1)) -» cf(r, r" | n c ). (36) 



4 2) (r,r" | an c ) ->■ 4 2) (r,r" | n c ), 



Then equations (27) and (28) become 



In (n c (r)A c 3 ) =/3 cy u c (r) + J dr'c® 



r, r' I n c )n c (r") 



\ng c {T,v' | n c ) = -/3 c c (r,r') + J dr"cf\r,r" | n c )n c (r") (g c (r",r' | n c ) - I) 



(37) 



(38) 



Together with the Ornstein-Zernicke equation (29) they are a closed set of equations to 



determine the density and pair correlation function for given potentials. This is the HNC of 
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liquid state theory, generalized to spatially inhomogeneous systems [21] . It is known to give 
very good results for uniform Coulomb systems [1J and for inhomogenous, confined Coulomb 
systems even at strong coupling conditions 



The corresponding inverse HNC forms for (30) and (31) 



PcfjL c (r) = In (n (r) A^) — J dr"c^(r, r" | n)n (r") , 



(3 c (p c (r, r') = — ln(l + h(r, r' \ n)) + h(r, r' | n) — c(r, r' | n c ), 
and the Ornstein-Zernicke equation (|3~lp,) is unchanged 

n) n fr") /i fr", r' I n 



c ( 2 ) ( r> r ' | n) = (r, r' | n) — J dr"c^ (r, r" 



(39) 
(40) 

(41) 



The "hole function", h(r,r' | n) has been introduced for notational simplicity 

h (r, r' \ n) — g (r, r' | n) — 1. 



(42) 



These equations, together with (34) for /3 c //3, provide practical forms to determine /3 c /j3, 
yS c // c (r), and ,9 c c (r,r / ) for the classical system, given appropriate quantum input. 



In the uniform limit (no external potential) the second term of (39) is simply related to 
the static structure factor S(k) 

1 



/WlnW/)../^-,-, 



n) 



(43) 



5(^ = 0)' 

For Coulomb systems, such as jellium considered in the following companion paper [19] and 
ideal Fermi fluids at zero temperature, S(k) vanishes as k — > and this contribution to /3 c /j, c 



diverges. Therefore, instead of (39) an alternative form obtained from a coupling constant 



integration and the HNC [25] is used, 



(3 c fi c = In (n c Aj?) -n J dr (c(r,n) + (3 c <p c (r) - ^h(r,n) (h(r,n) - c(r,n)) 



(44) 



III. THERMODYNAMICS - FURTHER CONSIDERATIONS 

The definition of the equivalent classical system assures that the pressure and density in 
the grand ensemble give the correct quantum results (in principle), e.g. 



PciPc I He, 4>c)V = -VL c (f3 c | /i c , <p c ) = p(/3 | fi, 4>)V = -Q(/3 | /i, 



(45) 
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However, some care is required in the definition of other thermodynamic properties. For the 
quantum system the thermodynamic variables are (3,/j,(r) for constant <p(r, r'). However, 
their variation leads to a change in all three classical variables /3 C , /x c (r), c (r, r'), since all 
are functions of /3 and functionals of /i(r). Therefore, the variation of the pressure leads to 
(in the following the volume V is always held constant) 



HpcV) 



dp c V 



d(3 c lMc,0c 



+ / drdr 



/ Sp c V , <90 c (r, r') 



+ / dr 



5 Pc V 
_S/i c (r) 



I /3c, 4>c 



+ / dr'dr 



cV r,r'j 9/3, 



y, r"l 



;r , r 



5fi c (r) 



5fx c (r) 



= S c dT c + J drn c (r)6fi c (r) (46) 

The second equality defines the classical thermodynamic entropy and thermodynamic den- 
sity in terms of the grand potential 



T r Sr 



P, 



dp c V 

d/3 c Uc ' 0c 



+ / drdr 



, Sp c V 
5<p c (r,r') 



I d(j) c (r, r') 



IMc' 



(47) 



n n r 



Sp c V 



+ / dr'dr' 



5p c V 



fr',r") 



> c (r',r") 



Similarly, the classical internal energy is defined by 

E c = T C S* C + / drn c (r)fjL c (r) - p c V 



5fi c (r) 



(48) 



dr/3 c /x c (r) 



5a 



5/3 c /j, c (r) 



df5 c 



l/?cMc 



(49) 



Note that derivatives in the last equalities of (50) and (49) do not have the restiction of 
constant <f) c . 

These same relationships hold for the quantum properties as well, since they are the 
general definitions of thermodynamics for the chosen variables, /J.,/3. In the quantum case 
is independent of the thermodynamic variables and hence is constant in the variations. 
This leads to the equivalent expressions in terms of equilibrium averages 

SQ , „ d(5VL . 



n(r) 



Mr) 



(n(r)> , E 



d/3 



H 



(50) 
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However, in the classical case the above leads to 



n c (r) = (n(r)) c - / rfr rfr — — |^ c (51) 



<50 c (r',r") IPc '^ <^ c (r) 

= + y " r " r k - k • (52) 

For example, the thermodynamic density, n c (r), differs from the average density of the text 
above, n c (r), because the latter is defined as a derivative of the grand potential at constant 

0c- 

IV. EXAMPLE - INHOMOGENEOUS IDEAL FERMI GAS 

To illustrate this definition of an equivalent classical statistical mechanics, consider the 
case of non-interacting Fermions in an external potential 

H-fiN K- J dr//(r)n(r). (53) 

Since H — jiN is now the sum of single particle operators, the pressure, density, and pair 
correlation function can be expressed in the form of single particle calculations 

p(f3 | n)V =p(/3\iJL,<j> = 0)V = f3'\2s + 1) J dr (r| In ^1 + e ~K^ (r) ) j | r ) , (54) 
n(r,/3|//)=n(r,/5|/x,0 = O) = (2s + l)(r| ^eK^~" (r) ) + |r) (55) 
g(r, r';/3 \ fi) = g(r, r';/3 \ fi, = 0) 

where (r| X |r') denotes a matrix element in coordinate representation, and s is the spin. 

Further reduction of these results to expose the dependence on (3 and \x (r) is difficult 
without solving the eigenvalue problem for the single particle Hamiltonian (p 2 /2m) — //(r). 
A useful practical approximation that captures most of the important exchange effects (but 
not any bound states if supported by the external potential) is obtained by replacing the 
operator /i(r) by its eigenvalue /x(r), a "local density approximation". The expectation 
values above then can be calculated as simple integrals 
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!/*)-»►— / rfr (2s + 1) 



h 3 



dp 



P_ 

2m 



-i 



n(r,/3 | /i) ->■ (2s + 1) h~ s / dp 



+ 1 



nr. r 



1 

h 3 



dpe 



rP'(r-r') 



e f(L-M(R)) + 1 



R 



(57) 
(58) 
(59) 



The results ( [571 ) are the familiar finite temperature Thomas-Fermi approximations for the 
thermodynamics, while (59) is its extension to structure [27J. The expressions for n(r,/3 | y) 



and n(r, r') are no longer functionals of /x(r), but rather local functions of /x(r) and /x(R) 



respectively. The change of variables from /3,/i to (3,n is accomplished by inverting (58) 
to obtain /i(r) = /i(r,/3 | n). Practical fits for n(r,/3 | //) and this inversion are recalled in 
Appendix [Bj along with simplification of n(r, r'). 

The corresponding classical results are non-trivial because (f) c (<li, <lj) ^ when 0(qj, qj) = 
0, because classical pair interactions are required to reproduce quantum exchange effects. 
Thus the thermodynamics and structure of a simple ideal quantum gas requires a corre- 
sponding classical system with the full complexity of an interacting many-body system. 
This classical many-body problem is addressed here and below in the HNC approximation 



described above by (39) and (40) together with the Ornstein - Zernicke equation (41) and 



(34) for f3 c //3. The solution for the classical parameters is a straight forward numerical task, 
but the results are different for each given external potential. To simplify the illustration 
here the following is further restricted to a uniform ideal Fermi gas (//(r) = //) in the ther- 
modynamic limit. See section 5 below for the non- uniform case of charges in a harmonic 
trap. For the uniform gas the equations for the classical system parameters simplify to 

1 







n 



6 



-n 



dvg{r)v ■ V/3 c c (r) 



P C H C = In (n c A^) - n j dr fc(r) + (3 t 



4>c{ r ) h(r) (h(r) — c(r)) 



p c </> c (r) = - In (1 + h(r)) + h(r) - c(r), 



c (r) = h (r) 



n 



J dr'c (|r — r'|) h (r') 



(60) 

(61) 
(62) 
(63) 



The superscript (2) on c^(r) and the dependence of c(r) ,h (r) on thermodynamic variables 
has been suppressed for simplicity. These last two equations can be solved for /3 c c (r) using 
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h(r) from (56) and (59) in the uniform limit. With that result, /3 c //3 can be calculated from 



(60), and then (3 c [i c determined from (44). Further elaboration is given in Appendix |B| 

The dimensionless potential /3 C C will be referred to as the Pauli potential. It is shown 
in Figure [l] as a function of the dimensionless coordinate r* = t/tq where r$ is the mean 
distance between particles defined by Autq/?) = 1/n. The state conditions are represented 
by r s = ro /as characterizing the density (where a# is the Bohr radius), and t = Pf/P for 
the temperature relative to the Fermi temperature (ftp 1 = ep = h 2 (37r 2 n) 2 ^ 3 /2m). For 
example, in these units n\ 3 = 8/ (37r 1/,2 t 3//2 ) and the classical limit occurs for t » 1 where 
the distance between particles is large compared to the thermal de Broglie wavelength. 
Figure [T] shows the Pauli potential at for t — 0, 1CT 1 , 1, and 10. 




FIG. 1: Ideal gas Pauli pair potential as a function of r* = r/r$ for t = 0, 0.1, 1, 10. 



Generally, the potential is positive, finite at r = (Pauli exclusion), and monotonically 
decreasing. The behavior is exponential at small r, but an r~ 2 algebraic tail develops for 



small t. This arises from the direct correlation function c^ 2 \r) in (62). Classical statistical 
mechanics does not exist for such a long range potential and it would appear that the 
equivalent classical system proposed here fails even for this simplest case of an ideal Fermi 
gas. However, this problem is "cured" for the corresponding case of classical Coulomb 
interactions with the same long range problem by adding a uniform neutralizing background 
(the one component plasma). The same procedure can be used here, i.e., a classical system is 
considered where the pair potential is supplemented in the Hamiltonian by a corresponding 



uniform compensating background. The pressure equation (26) is modified due to this 
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background by the replacement of g c (r,r') by g c (r, r') — 1, and (60) becomes 



f3c _ n_ 



1 n / drh(r)r ■ V/3 c c (r) 

6 



(64) 




FIG. 2: Ideal gas reduced classical temperature t c i = T c i/Tp as a function of t = T/Tp ■ Also 
shown is the result of PDW. 

Figure |2] shows the classical temperature relative to the Fermi temperature, f3p/(3 c = t c , 



as a function of t obtained from (64). It is seen that the classical temperature T c remains 
finite at T = in all cases, and crosses over to T c = T at high temperatures. The PDW 

1 /2 

model postulates the form T c = (T 2 + TJ 2 ) . The model originally uses the average energy 



per particle at T = to evaluate T = 2Tp/h. The result from (64) is quite close T c (t — 0) ~ 
0.43. To compare the dependence at finite t, the PDW form is also shown in Figure [2] with 
T = T c (t — 0). It is seen that the results are quite similar although the PDW form has a 
somewhat faster cross over to the classical limit. 



Figure [3] shows similar results for \x c j Ep as a function of t obtained from (61 ). Also shown 
is the result for the quantum fi/Ep. Both forms depend only on t (independent of r s ). At 
high temperatures the chemical potential of the representative system goes over to ln(nA 3 ) 



as in equation (61 ) while the quantum chemical potential goes over to ln(nA /2). They differ 
by a factor of ln(2) at high temperatures due to the internal spin degrees of freedom in the 
quantum calculation. This is accommodated in the comparison shown in Figure |3j 

With the parameters of the equivalent classical system determined, it is now possible 
to calculate any desired property of the ideal quantum gas by classical methods. Since 
the pressure is given exactly by the definition of the classical - quantum correspondence 
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-4- 

-5 



-H/E p +tln(2) 
-ME, 



t=T/T 



FIG. 3: Ideal gas dimensionless chemical potential /J, c /Ep as a function of t. Also shown is the 
corresponding quantum chemical potential jx/Ep + iln2. 

equations, it is useful to consider the internal energy. First, note that all the dimensionless 
thermodynamic quantities of both the quantum and classical ideal gas depend on only one 
thermodynamic parameter z = e^ M . In particular 

PcPX = G(z). (65) 

The specific form for G(z) is not required for the present discussion. The corresponding 
quantum result is /3 c p c \l = 73/2(2), where 73/2(2) is the Fermi integral of Appendix [BJ 



(B19). Although G(z) and 73/2(2) are quite different, the fact that they both depend only 
on z implies that the relationship among different thermodynamic properties is the same in 
both classical and quantum cases. For example, the energy per particle is determined from 



the pressure via the thermodynamic definition (49) with (61) 



~ <9A -3 2 

E d = -VG(z)^-\ z = 3 p d V (66) 

The contribution from dG(z)/d(3 c \ Zc y vanishes since z = z(z c i). This follows on dimen- 
sional grounds since z is dimensionless, and there is no additional energy scale to make /3 C ; 
dimensionless. Hence the classical calculation gives the known exact quantum result. This 
non-trivial result for a classical interacting system is a strong confirmation of the effective 
classical map defined here. As noted in the last section, it can be verified that the classical 
average of the classical Hamiltonian does not give the correct relationship to the pressure. 



Instead it is necessary to define the internal energy thermodynamically, as is done in (66) 

19 



In a similar way it is verified that the relationship of the classical entropy to the pressure 

and density is the same as in the quantum case 

5 

T d S d = -p d V - fi d N d . (67) 

This is a form of the Sackur- Tetrode equation valid for the quantum case. It is emphasized 
that these simple ideal quantum gas results are being retained for the more complex effective 
classical system with pair interactions via the Pauli potential. 



V. WEAK COUPLING PAIR POTENTIAL 



For systems with real interactions 0(qi,q2) between the particles, the classical pair 
potential will have the form 

&0c(qi,qa) = (<k(qi, q 2 )) (0) + A(q 1; q 2 ), (68) 

where (0 c (q 1; q 2 ))^ ^ is the ideal gas Pauli potential and A(q 1 ,q 2 ) denotes the contribution 
to the effective potential from the real pair potential. Obviously, in the classical limit 
A(qi, q 2 ) — > /30(qi, q 2 ). Another exact limit is the weak coupling limit for which the direct 
correlation function becomes proportional to the potential pQ, or stated inversely, 

PMr) -> -c(r), (/3 c c (r)) (o) -> -c^(r). (69) 
Thus a possible approximation incorporating this limit is 

PMr) -> (/3 c c (r)) (o) - (c(r) - c^(r)) (w) , (70) 
where (c(r) — c^(r))^ w denotes a weak coupling calculation of the direct correlation func- 



tions from the Ornstein - Zernicke equation (29) 



(/3 c c (r)) W = (r)+n dv' (/3 c c (|r - v'\)) {w) h^(r'), (71) 



with h w (r) being the quantum hole function in its weak coupling form. In the following 
companion paper this approximation is applied to the electron gas (jellium) for which the 
weak coupling form is given by the random phase approximation. 



To further interpret approximation (70) consider a uniform system and rewrite 



A(qi, q 2 ) = A(|qi — q 2 |) in terms of the static structure factor 

S{k) = 1 + n J rfre ik - r /i(r), (72) 
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= M\<il ~ Qal) " / rfqG(q)/30(| qi - q 2 - q|) (73) 

Here G(q) is known as the local field corrections in linear response for the classical system. 
It is seen that the quantum corrections to |qi — q 2 1 ) can be interpreted as local field 
corrections [20, EE], although due to quantum effects rather than classical correlations. In 



the case of jellium, the qualitative features of (73) are a regularization of the Coulomb 
singularity at qi = q 2 , and cross-over to an asymptotic Coulomb decay for large |qi — q 2 | 
with effective coupling constant given by the exact perfect screening sum rule. 



VI. DISCUSSION 



The objective here has been to define an effective classical equilibrium statistical mechan- 
ics that corresponds to a chosen quantum system of interest. The motivation is to allow 
application of existing strong coupling classical methods (e.g. liquid state theory, MD sim- 
ulation) to calculate properties of the quantum system under conditions for which current 
theoretical approaches are not adequate. The correspondence has been defined by equiva- 
lence of the pressures, densities, and pair correlation functions for the classical and quantum 
systems. In this way, the relevant parameters for the classical grand ensemble are fixed - 
the temperature, local chemical potential, and pair potential. These classical parameters 
were given a more explicit representation in terms of the quantum parameters by inverting 
the classical many-body problem using the HNC integral equations. Formal questions such 
as the existence of this inversion have not been addressed, and only presumed to hold. A 
counter example is given by the classical representation of jellium, where the equivalence of 
pressures is not possible under conditions of negative pressures. 

The three correspondence conditions of section [II] are not unique, and other choices may 
be preferred in specific applications. Furthermore, applications require the introduction of 
appropriate approximate forms for these correspondence conditions that should be tailored 
to the particular system at hand. The special case of a uniform ideal Fermi gas was illustrated 
using the HNC integral equation to determine the parameters of the effective classical system. 
A peculiarity is the long range nature of the effective pair potential at very low temperatures, 
requiring the introduction of a compensating uniform background. The resulting classical 
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thermodynamics was shown to reproduce the exact relationships of various thermodynamic 
functionals (e.g., pressure, internal energy, entropy). 

The non-uniform ideal Fermi gas is more interesting and its representation as a func- 
tional of the local density is a fundamental problem within density functional theory [24|. 
Approximations such as the Thomas-Fermi representation have limited applicability. An 
effective classical representation along the lines described here would provide access to bet- 
ter approximations, since the functional dependence on density for the classical system is 
simple. 

A weak coupling approximation for the effective pair potential in systems with real forces 
was described in the last section. With that pair potential known, the effective chemical 
potential and effective temperature can be calculated. This is illustrated in the following 
companion paper for the uniform electron gas [19] . The pair correlation function is calculated 
for a wide range of densities and temperatures, and good agreement is obtained with diffusion 
Monte Carlo results at zero temperature. 

A second application in that paper is to shell structure for charges confined in a harmonic 
trap. Classically, shell structure arises only from Coulomb correlations [22J. A preliminary 
investigation there shows new origins of shell structure due to diffraction and exchange, even 
in the absence of Coulomb correlations (mean field approximation). 



VII. ACKNOWLEDGEMENTS 



This research has been supported by NSF/DOE Partnership in Basic Plasma Science and 
Engineering award DE-FG02-07ER54946 and by US DOE Grant |de-sc0002139l 



Appendix A: Exact Coupled Equations for n c (r) and g c (r, r') 



The objective of this appendix is to outline the origin of the exact equations (27) - (29) for 
the classical density and pair correlation function. First, make a change of thermodynamic 
variables from temperature and chemical potential /3 c ,/i c (r) to temperature and density 
(3 C , n c (r) by the Legendre transformation 

F C {/3 C | n c ) = tt c (f3 c | /i c ) + J dvfi c (r) n c (r) , (Al) 



22 



where now the free energy is a functional of the classical density rather than \x c (r). Their 
relationship is given by the first derivative 



SF C 
5nJr) 



(A2) 



Here, and throughout this Appendix derivatives are taken at constant <ft c , so the densities 



involved are those defined as in (17). The free energy is now divided into its ideal gas 
contribution = —f3^ 1 J dr [1 — In (n c (r) Aj!)] n c (r), where A c = {2nh 2 f3 c /m) 1 ^ 2 , and the 



remainder F CfiX (excess free energy), so that (A2) becomes an equation for the density 



In (n c (r; C ) \ 3 C ) = /3 c /i c (r) + c« (r; e \ n c ) . 



(A3) 



Here c£ (r | n c ) is the first of a family of functions (direct correlation functions) defined by 
derivatives of the free energy 

5 m F 



ri, ..,r m ;P c I n c ) = -(3 



Sn c (r 1 )..Sn c (r r 



(A4) 



Equation (A3) relates the density n c (r) to a given external potential (recall /i c (r) = 
A*c — 4>c,ext,( r )) ■ Consider now a different external potential given by 4> c ,ext,( r ) + 0c( r > r ') and 
associated density n c (r, r'). This external potential corresponds to the original one 4> c ,ext,{ r ) 
plus a new source of potential of the same form as would occur if another particle were 
added at the point r'. It follows that this new density is proportional to the pair correlation 
function for the original system pQ 



n c (r, r) = n c (r) g c (r, r') . 



(A5) 



The equation corresponding to (A3) for this new external potential is 



In (n c (r) g c (r, r') A^) = /3 c /i c (r) - (3 c (fi c (r, r) + (r; (3 C \ n c g c 



(A6) 



Finally, subtracting (A3) from (A6) gives the desired equation for g c (r, r') 



In (g c (r, r') A*) = -/3 c c (r, r') + c« (r; f3 c \ n c g c ) - c« (r; (3 C \ n c ) 



(A7) 



The notation used implies that the functional (r; /3 C \ •) in both (A3) and (A6) are the 
same. This follows from density functional theory where it is demonstrated that the free 
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energy is a universal functional of the density, the same for all external potentials. Equations 



(27) and (28) now follow directly from (A3) and (A7) and the identity 



c« (r; p c \X) = C W (r; {3 C \Y)+ [ dad a c^ (r; c \aX + (l- a) Y) 

Jo 

,5c ( c } (r;& \aX + (l-a) Y) 



c« ( r; P B \Y)+ I da I dr 



5n r (V 



(X (r') - Y (r')) 



4 X) (r; C \Y) + j da j dr'cf (r, r'; (3 C \ aX + (1 - a) Y) {X (r') - F (r)) 

(A8) 



with appropriate choices for X and Y. 



The Ornstein - Zernicke equation (29) is an identity obtained as follows. The second 



functional derivative of the grand potential is related to the pair correlation function by 



5 2 (-(3M 6n c (r>) 



5fi c (r) (i c (r') 5(i e (r) 
Similarly, the second derivative of the free energy is 



f3 c n c (r') [S (r - r') + n c (r) (g c (r, r') - 1)] , (A9) 



5 2 F r . 



5n c (r') 5n c (r) 5n c (r') 
Then the chain rule 



can be written 



S = A"V(')-[«('-r') 



n r (r) cl 2 -* (r, r' | n c 



dr' 



Sn c (r) S/j, c (r") 
8n c (r") <5n c (r') 



5 (r - r'l 



(A10) 



(AH) 



J dr" [5 (r - r") + n c (r) (g c (r, r") - 1)] [5 (r"—r') - n c (r") 4 2) (r", r' 



ra c )] = 5 (r - r') . 

(A12) 



This gives the Ornstein-Zernicke equation (29). 



Appendix B: Inhomogeneous Ideal Fermi Gas 

The thermodynamic and structural properties of an inhomogeneous ideal Fermi gas are 
straightforward to calculate in a representation that diagonalizes the effective single particle 
Hamiltonian 

|^-M?))^k(r) = e k ^ k (r), (Bl) 
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where k labels the corresponding quantum numbers. For Fermions with spin s, the quantum 
numbers are labeled by k — (s, k). The Hamiltonian in second quantized form is then simply 

H = J2^{a K , (B2) 

where a) K) a K are the creation and annihilation operators for occupation of the states {ipk}. 
Then the pressure is found directly from evaluation of the grand potential Q c 



p(/3 | n)V = /T 1 + e_/3ek ) = P~ l Tr {l) In 1 + e 

K \ 

= (2s + l)r i y dr(r|ln ^1 + e"K^" M?) ) j | r ) , 



(B3) 



where a coordinate representation has been used in the last expression. 

The local density and pair correlation function are obtained from the one an two particle 
density matricies. In the diagonal representation these are 

P (1) («i; « 2 ) = (4i a K2> = (4i a /«i) <W 2 ( B4 ) 

P {2) («!, K 2 ; « 3 , K4) = (4i a lt, a «4 a «s) = ( 5 «1A3 5 K2A4 - ^1,«4 5 «2,« 3 ) (4i a «l) (4 2 a K 2 ) ( B5 ) 

Where the mean occupation number is 

(a{a K ,) = 5 K , K ,(e^ + l)~\ (B6) 
The coordinate representations are 

pW (r,(7 l5 r', a 2 ) = £ ^ (r) ^ (r') «<»„) = 5 n „ (r| In (V^"^) + l") |r'> 

ki,k 2 V / 

= <5 CTl , CT2 n(r,r') (B7) 



p (2) (n, (Ji, r 2 , a 2 ; r' 1; a 3 , r' 2 , a 4 ) = ^ ^ ( r i) ( r 2) Vv 3 ( r i) ( r 2 ) ^X^^a} 

ki..k6 

= 8„ U V3 S *2,V4 ( ri ) ^ ( r 'i) (4 1 °«i> $^C 2 ( r 2)^« 2 (r 2 ) <4 2 o K2 > 

ki k 2 

- <W 4 <W 3 X^«i ( ri )^«i ( r 2) (4i a «i>5^^«a ( r 2) ^ K2 (ri) (a{ 2 a K2 ) 

ki k 2 

(B8) 



The diagonal elements are 



p (1) (r,<ri;r,<7i) = n(r,r) 



(B9) 



p (2) (n, oi, r 2 , cr 2 ; n, 01, r 2 , a 2 ) = p (1) (ri,<7i; r 1; a x ) p (1) (r 2 ,a 2 ; r 2 , a 2 ) 



<W 2 P (1) (ri,cn; r 2 , cri) p (1) (r 2 ,a 2 ; r 1( er 2 ) (BIO) 



Finally, the density and pair correlation function are identified from the summation over 
spin states 

n(r) = P (1) (r^U r, 0i) = (2s + 1) n (r, r) (Bll) 



n 



ri)n(r 2 )^ (ri, r 2 ) = p (2) (n, oi, r 2 , cr 2 ; n, o"i, r 2 , a 2 ) = n(ri)n(r 2 )-(2s + 1) n (r x , r 2 ) n (r 2 , n) 

(B12) 



This gives the results (55) and (56). 



The local density and pair correlation function are determined from the function n (r, r') 



obtained from the single particle density matrix (B7) 



n(r, r') = (r| ( e 



1 10 



(B13) 



In the local density approximation of the text, p(r) — > /u(R), where R = (r + r') /2, this 

(B14) 



becomes (59) 



Further simplification is possible to get 

2A ro ° 

with 



n(r,r')A 3 = — -/ tfaa: ( ^(Rje* 2 + 1 ) sin (2>/nx |r - r'l /A) 

7r r — r l/ V / v 



*(R) = e^ (R) 
Accordingly the density and pressure simplify to 



n(r)A 3 = (2s + 1) n(r, r)A 3 = (2s + 1) / 3/2 (*(r)), 



/3pA a = - / dr(2 S + l)/ 5/2 (z(r)), 



/a/ 2 (^) = ->= / dxx 2 (z 1 e x +1) , / 5/2 (;z) = 



(B15) 
(B16) 

(B17) 
(B18) 



with the definitions 

4 r°° 

'o 



cfex 4 [2; V + 1 



(B19) 
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